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The problem of two-dimensional steady nonlinear dynamics in plane Couette flow is revisited 
using homotopy from either plane Poiseuille flow or from plane Couette flow perturbed by a small 
symmetry-preserving identity operator. Our results show that it is not possible to obtain the non- 
linear plane Couette flow solutions reported by Cherhabili and Ehrenstein [Eur. J. Mech. B/Fluids, 
14, 667 (1995)] using their Poiseuille-Couette homotopy. We also demonstrate that the steady solu- 
tions obtained by Mehta and Healey [Phys. Fluids, 17, 4108 (2005)] for small symmetry-preserving 
perturbations are influenced by an artefact of the modifled system of equations used in their paper. 
However, using a modifled version of their model does not help to flnd plane Couette flow solution 
in the limit of vanishing symmetry-preserving perturbations either. The issue of the existence of 
two-dimensional nonlinear steady states in plane Couette flow remains unsettled. 

PACS numbers: 47.20.Ft, 47.15.Fe, 47.20.Ky, 47.27.Cn, 47.27.De, 47.27.N-, 47.32.Ef 



I. INTRODUCTION 

Even though plane Couette flow is known to be lin- 
early stable for all values of the Reynolds number Re, 
transition to turbulence is observed experimentally and 
numerically at quite modest values of Re ~ 350 in this 
flowi. As a consequence, the quest for nonhnear instabil- 
ities has been the subject of many studies. Since there is 
no direct connection between the basic state and poten- 
tial nonlinear branches of solutions, homotopy (continu- 
ation) from other flows subject to Hnear instabilities has 
been a very popular method to discover such branches. 
Exact three-dimensional solutions have been reported by 
Nagata^ using homotopy from linearly unstable rotating 
plane Couette flow, by Clever and Bussed using homo- 
topy from a convectively unstable plane Couette flow and 
by Waleffe^'— using a forcing homotopy. These three- 
dimensional states match very well experimental data on 
transition to turbulence in plane Couette flow, but other 
solutions, including two-dimensional ones, may well be 
participating in phase-space dynamics and in the transi- 
tion process. 

Another kind of homotopy was suggested by Mili- 
nazzo and Saffman— to discover two-dimensional nonlin- 
ear solutions in plane Couette flow by connecting the 
nonbifurcating no-slip plane Couette flow to the no-slip 
plane Poiseuille flow, which becomes linearly unstable at 
Re = 5772 (Zahn et al.-i, Herbert^). Starting from the 
traveUing-wave solutions in this flow and progressively 
changing the original background flow from a Poiseuille 
to a Couette proflle, Milinazzo and Saffman— attempted 
to continue these solutions to a pure Couette flow, but 
they failed to achieve this goal. Their lack of success may 
have to do with resolution: Cherhabili and Ehrenstein— 
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(hereinafter CE95) used the same approach with more 
Fourier modes in the streamwise direction and claimed to 
have identifled steady two-dimensional locaHzed nonhn- 
ear solutions for the Couette flow. Cherhabili and Ehren- 
stein— subsequently computed the stability of these two- 
dimensional steady states and reported the discovery of 
three-dimensional nonlinear structures which are differ- 
ent from those found by Nagatai^. 

Recently, Mehta and Healey— (hereinafter MH05) sug- 
gested yet another interesting homotopy to compute two- 
dimensional nonlinear states in plane Couette flow. They 
argued that the homotopy from Poiseuille to Couette 
flow was not symmetry-preserving, so that possible non- 
linear states relying on the 0{2) symmetry of the Cou- 
ette flow might not be found by continuation from flows 
that do not have this symmetry. Instead, they proposed 
a symmetry-preserving homotopy based on small per- 
turbations of the linearization of the Couette operator 
that causes the least stable eigenvalues of the plane Cou- 
ette flow to become unstable. Doing so, they hoped to 
be able to track the nonlinear branch associated with 
this bifurcation point down to zero perturbation and to 
obtain two-dimensional nonlinear states for the unper- 
turbed plane Couette flow. For small perturbation pa- 
rameters, they discovered steady solutions looking similar 
to the solutions of CE95. MH05 notably describe the nu- 
merical appearance of so-called "frozen- waves" with phase 
speed c = in perturbed Couette flow starting from bi- 
furcating traveUing-wave solutions at Re > 390. CE95 
had also noticed the numerical appearance of solutions 
with c = in their calculations. However, these solu- 
tions cannot be followed down to zero perturbation (ex- 
act plane Couette flow) , which raises the questions of the 
existence of the solutions reported by CE95 and of the 
connections between their flndings and those of MHOS. 

Considering that the case is quite unclear, we attempt 
to address some of these issues in this study. We flrst 
revisit the Poiseuille-Couette homotopy problem and ar- 
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gue that the solutions obtained by CE95 are spurious. 
Then, we clarify the nature of the steady solutions found 
by MHOS at large Re for small values of their pertur- 
bation parameter e. We show that a direct connection 
exists between travelling and steady solutions for the per- 
turbed flow and demonstrate that their high Reynolds 
number steady solutions are influenced by an artifact in 
the model that introduces a mean flow eigenmode. We 
obtain solutions for small perturbations using a revised 
symmetry-preserving homotopy and show that they can- 
not be continued to non-perturbed plane Couette flow 
either. Our results therefore suggest that it is not pos- 
sible to identify two-dimensional nonlinear solutions in 
plane Couette flow using these two particular methods. 

The paper is organized as folllows. In Sec. Ull we 
formulate the physical problem, describe the underly- 
ing symmetries and present the numerical techniques 
used to compute nonlinear branches of solutions. Sec- 
tions [ini and [IV] describe our efforts to obtain steady 
two-dimensional nonlinear plane Couette solutions using 
either the Poiseuille-Couette homotopy with the same 
procedure as CE95 or the symmetry-preserving homo- 
topy of MH05. A short discussion of the results (Sec. |V| 
concludes the paper. 



the nondimensional equations governing the evolution of 
U and P are the incompressible Navier-Stokes equations 

dtU + U-VU = -VP+^AU , V-U = Q. (2) 

Rc 

They are supplemented by periodic boundary conditions 
in X and by no-slip boundary conditions in y for u and v: 

u{x, — 1, = u{x, 1, t) = v(x, —1, t) = l'(.T, 1, = , 

(3) 

For numerical purposes, a Fourier decomposition is used 
in the x direction and a scalar equation for the shear- 
wise velocity v (or equivalently for a stream function V') 
governing the evolution of the ^ Q modes is solved: 

(^dt - Ai; + • ((7 • V(7) - , (4) 

where P„ = — ej,-VxVx(-). An auxiliary equation 
describing the evolution of the mean streamwise velocity 
fleld must be solved simultaneously for completeness: 

Rc {dtu + dymi) — dyll , (5) 

where overbars denote x-averaged quantities. 



II. EQUATIONS, SYMMETRIES, NUMERICS 
A. Physical model and background state 

In the following, an incompressible fluid with kinematic 
viscosity v and constant density p taken to be unity 
is considered. The plane shear flows investigated here 
are independent of the streamwise coordinate x and are 
bounded by two rigid walls in the shearwise y direction, 
with —d < y < d, where d is chosen as length unit. The 
X direction is taken to be periodic with spatial period 
Lx = 2Trd/a. A reference velocity Uo is deflned and a 
mean streamwise pressure gradient —21^ Uo{l''ii)/d^ and 
mean streamwise wall velocities Usiy = ±1) = ±?7 Uo are 
imposed in order to maintain the following generic (non- 
dimensional) streamwise oriented background shear flow: 

UB{y) = ii-v){i^y^)+vy ■ (i) 

Here, 77 is a nondimensional parameter that allows contin- 
uation from plane Poiseuille flow (77 = 0) to plane Couette 
flow (r/ = 1). The Reynolds number of the flow, based on 
Uo and half of the channel width, reads Re = {Uod)/v. 



B. Equations 

The full non-dimensional velocity fleld, denoted by U , 
is decomposed into the original streamwise background 
flow (C/s,0) and two-dimensional velocity perturbations 
denoted by m = (w, v). The corresponding total pressure 
is P = Pb + P, where p is a pressure perturbation, and 



C. Symmetries 

Depending on the value of r/, this set of equations has 
different symmetries. For ?/ 7^ 1, the only symmetry is 
the translation symmetry 5*0(2) 

SO{2) : {x,y) ^ {x + a,y) , {u,v,p) ^ {u, v,p) (6) 

with a e [0, 27r/a]. For plane Couette flow (77 ~ 1), an 
extra (point) reflection symmetry is present, 

: (x,?/) (-X, -y) , {u,v,p) {'~u,-v,p) . (7) 

In that case, the full flow has the more general 0{2) sym- 
metry. 

D. A symmetry preserving homotopy 

An extra symmetry-preserving artiflcial term propor- 
tional to eu, where e is a small non-dimensional parame- 
ter, can be added to the r.h.s. of the original Navier- 
Stokes equations ^ in order to cause the least sta- 
ble eigenvalue to become unstable and to perform the 
symmetry-preserving homotopy mentioned in the intro- 
duction. This technique is described in detail in MH05. 

E. Numerics 

The exact two-dimensional nonlinear steady and trav- 
elling solutions of Eqs. ©-([S]) presented in the next sec- 
tions have been obtained thanks to a nonhnear continu- 
ation codoi^ similar to that of Waleffe^ and Wedin and 
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Kerswell— . Chebyshev collocation on a Gauss-Lobatto 
grid in the shearwise direction and a Fourier represen- 
tation in the x direction are used. Clamped boundary 
conditions are implemented using differentiation matri- 
ces computed with the DMSuite packagei^. Nonlinear 
terms are computed in physical space and steady non- 
linear equations are solved in spectral-physical-spectral 
space using Newton iteration. A direct LAPACK 
solver is called at each Newton iteration and dealias- 
ing is performed in the a;-direction. Pseudo-arclength 
continuatioiii^ is implemented to follow branches of solu- 
tions. The code has been tested quantitatively by solving 
the classical problem of subcritical continuation of two- 
dimensional traveUing waves in no-sHp plane Poiseuille 
fiow^i^. The full three-dimensional version (not used 
here) also quantitatively reproduces the continuation di- 
agrams of three-dimensional nonlinear solutions obtained 
by Waleffe^ in no-slip, non-rotating plane Couette fiow^^. 



III. POISEUILLE-COUETTE HOMOTOPY 

The Poiseuille-Couette homotopy was first considered 
in order to compute the two-dimensional nonlinear so- 
lutions reported by CE95 in plane Couette flow. No 
particular assumption regarding the symmetries of the 
solutions was made in order to perform these compu- 
tations and the numerical method employed here copes 
with solutions with either zero or non-zero phase speed 
(the transformation dt = —cdx, where c is the phase 
speed, is used in the latter case). All our attempts proved 
unsuccessful. We simply followed their procedure and 
performed continuation with respect to 77, starting from 
two-dimensional nonlinear travelling wave solutions in 
Poiseuille flow (ry = 0) at Re = 15993, which corresponds 
to the flrst case considered by CE95. Continuation curves 
for mixed Poiseuille-Couette flow are shown in Fig.[TJ For 
the maximum affordable resolution and a = 1 , a turning 
point is found at 77 ~ 0.23, therefore it proves impos- 
sible to isolate solutions for plane Couette flow. Since 
CE95 found Couette solutions for small a only (albeit 
for smaller Re) , we tried to perform further continuation 
with respect to a for various values of 77 and different res- 
olutions including the maximum resolution used by CE95 
(corresponding to = 52 in our code). The results, pre- 
sented in Fig. m show that it is not possible to identify 
plane Couette flow solutions for small values of a either. 
The nose of the surface of solutions in the a — i]— wall 
shear rate space is located at 77 ~ 0.5 and a ~ 0.2 for 
Rc — 15993, even for rather large streamwise resolution. 
This result is similar to Fig. 1 of CE95, who presented 
low streamwise resolution solutions down to 77 = 0.45 for 
Re = 15993. As already acknowledged in CE95, it is clear 
from the flgure that there are convergence issues at low 
a for such a high value of Re, even with a resolution such 
as {Nx,Ny) = (52,48). We notably encountered difficul- 
ties performing the low-a continuation at 7; = 0.45. The 
same procedure was consequently attempted at the less 
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FIG. 1: Continuation of two-dimensional solutions with a = 1 
in mixed Poiseuille-Couette flow (HJ as a function of 77, for 
Re = 15993 (black curves with turning points at 77 ~ 0.23 
and Re — 7003 (blue online, with turning points at 7? ~ 0.13). 
Full line: {Nx,Ny) — (52,48), dashed line (almost undistin- 
guishable from full line): {N:^,Ny) = (32,48), dotted line: 
{Nx,Ny) = (6,48). The non-dimensional wall shear rate 
\dU /dy\(y = ±1) has been used as a measure of nonlinearity 
(for the laminar flow, \dU/dy\{±l) = 1). The "gap" at ?? ~ 
results from the existence of a turning point close to 77 ~ 0.02. 



resolution-demanding Re = 7003 at which CE95 found 
their Couette fiow solutions. The corresponding contin- 
uation curves for mixed Poiseuille-Couette fiow are also 
shown in Fig. [TJ This time, a turning point is found at 
77 ~ 0.14 for Q = 1. Further continuations with respect 
to a for various values of rj and the same maximum reso- 
lution as the one used by CE95 are represented in Fig. [3] 
and show that it is not possible to identify plane Couette 
fiow solutions for small values of a either. The nose of the 
surface of solutions in the a — 77— wall shear rate space is 
clearly located at 7; ~ 0.14 and a ~ 1 for Re = 7003. The 
same procedure was finally attempted for even smaller 
values of Rc and led to similar negative results. 

The main trend is therefore that it seems possible to 
reach solutions with lower 77 with increasing Re, but it is 
clear from our results that no 77 = 1 solution can be found 
for the range of Re investigated by CE95. Since only very 
little information is given on their precise procedure to 
obtain two-dimensional solutions, it is actually very hard 
to figure out where the problem could come from in their 
results. We originally thought that aliasing errors may 
be present in their calculations, but since the numeri- 
cal method used by CE95 is fully spectral, their results 
should not in principle suffer from this problem, at least 
in the streamwise direction. One possibility would be 
that their error tolerance for the Newton algorithm was 
not small enough. In the present calculations, conver- 
gence was considered to be achieved for an error smaller 
than 10~^ in energy and a precision of 10~^ was reached 
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FIG. 2: Continuation of two-dimensional solutions (Re = 
15993) in mixed Poiseuille-Couette flow ([T|) as a function of a 
for 77 = 0.17 (black curves with the largest a values), r] = 0.23 
(blue online, curves centered on a ~ 0.5), r] = 0.45 (red online, 
narrow curves at low a. Full line: (A''^,, A^^) = (52, 48), dashed 
line: (iV^.iVy) = (32,48), dotted \ine:{N^, Ny) = (6,48). So- 
lutions for A'^j; = 6 do not extend to 77 = 0.45 and the corre- 
sponding red dotted curve at a ~ 0.2 is for 77 = 0.4. Note the 
convergence problems at low a for {Nx,Ny) = (32,48) and 
(iV,,iV,) = (52,48). 



in many cases. To conclude this section, it is worth not- 
ing that the Milinazzo and Saffman— and CE95 homo- 
topy ^ is only a particular choice of Poiseuille-Couette 
homotopy, and that more general transformations from 
one flow to the other, e. g. 

UB{y)^{i-.f{v)){i-y^)+9{v)y , (8) 

where / and g satisfy the right end conditions, may lead 
to different results. This kind of behaviour has notably 
been reported in Taylor-Couette flow by Rucklidge and 
Champneys^, who discovered that a connection between 
solutions for two different sets of boundary conditions 
was only possible in their problem for some special ho- 
motopy choices. 



IV. SYMMETRY-PRESERVING HOMOTOPY 

Motivated by the apparent similarities between the 
steady solutions of CE95 and those that MH05 ob- 
tained for small values of their perturbation parameter 
e, we then considered the problem of finding steady two- 
dimensional nonlinear solutions in plane Couette fiow 
{i] ~ 1) using the homotopy proposed by MH05 and 
parameters similar to those used in their calculations. 
An s symmetry-preserving identity perturbation was in- 
cluded in the equations for that purpose. As mentioned 
in Sec. m this causes the least stable eigenvalue of the 



FIG. 3: Continuation of two-dimensional solutions (Re = 
7003) in mixed Poiseuille-Couette flow ((l]) as a function of 
a for g = 0.05 (black, family of large curves), 77 = 0.12 (blue 
online, family of medium size curves), 77 = 0.132 (red on- 
line, family of small curves corresponding to the "nose" of the 
surface of solutions). Full line: {Nx,Ny) — (52,48), dashed 
line (undistinguishable from full line curves): {Nx,Ny) = 
(32,48), dotted \me:{N^,Ny) = (6,48). The results for 
{Nx,Ny) = (52,48) presented here are also undistinguishable 
from those (not shown here) obtained for the same resolution 
{N:,,Ny) = (52, 28) as that used by CE95. 



Couette fiow to become unstable. For Re > 375 and 
a = 0.2, the bifurcation is of the Hopf type, while for 
Rc < 375 and the same a, the bifurcation is steady. In 
the latter case, for an identity perturbation with e > eo, 
where eo — 0.0456, the eigenvalue crosses the imaginary 
axis for two values of Re (see MH05), leading to two 
different branches of steady nonlinear solutions (Fig.IH. 
MHOS showed that steady branches with e < eq can- 
not be computed directly because they reduce to a single 
point at e = eo and Re ^ 270 (see their Fig. 11). Never- 
theless, as noted in their paper, there is a connection be- 
tween low Re branches with e > eo and high Rc branches, 
including some with e < eo. For instance, a solution at 
e = 0.04 (dashed line) can be obtained by reducing e, 
starting from a high energy solution at e = 0.05 and 
Re ~ 360. The associated e = 0.04 steady branch can 
then be continued to large values of Rc and transformed 
into high Re branches with different values of e. 

In Fig. m a travelling branch at e = 0.04547 is also 
plotted. This branch connects to a steady high Rc branch 
with the same e at Re = 500. The presence of such a trav- 
elling branch, which was discovered by performing con- 
tinuation with respect to e in a range of Rc where steady 
branches of solutions also exist, confirms the findings of 
MH05 that travelling solutions bifurcating at Re > 375 
can turn into steady solutions by changing e at fixed Re. 
Fig. [5] shows a typical continuation with respect to e that 
clarifies this connection. Continuation in this figure is 
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FIG. 4: Bifurcations and nonlinear continuation of two- 
dimensional solutions of the £-perturbed plane Couette flow 
of MH05 with a = 0.2, e = 0.05 (full line), e = 0.04 
(dashed line) and e = 0.04547 (dotted and dash-dotted lines). 
{Nx,Ny) = (32,48) has been used (the results of some of these 
computations done at (Nx,Ny) = (52,48) are undistinguish- 
able from those at A^'^, — 32). All solutions are steady except 
for the dash-dotted line branch at e = 0.04547 (dash-dotted, 
red online), for which the phase speed reaches when it meets 
the dotted steady branches. The vertical dash-triple dotted 
line corresponds to the path followed in the continuation pre- 
sented in Fig. [S] Continuation was stopped at large energy 
(strong nonlinearity) due to insufficient numerical resolution. 
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FIG. 5: Nonlinear continuation of two-dimensional solutions 
of the £-perturbed plane Couette flow of MH05 with respect 
to e, for a = 0.2, Re = 500 and {Nx,Ny) = (32,48). The 
dashed line corresponds to the phase speed of the solutions 
and the full line to their kinetic energy. The starting point 
of the continuation corresponds to the Hopf bifurcation of 
the perturbed Couette flow, which occurs at e = 0.05358 for 
Re = 500. Continuation was stopped at large energy (strong 
nonlinearity) due to insufficient numerical resolution. 



FIG. 6: Bifurcations and nonlinear continuation of steady 
two-dimensional a = 0.2 solutions of the e-perturbed plane 
Couette flow without the artifact mean flow perturbation for 
e = 0.05 (full line), e = 0.05358 (dashed line) and e = 0.04547 
(dash-dotted line) when the e perturbation term is added to 
the fci / part of the equations only. {Nx,Ny) = (32,48). 
Continuation was stopped at large energy (strong nonlinear- 
ity) due to insufficient numerical resolution. 



done along the line Re = 500 represented in Fig. |4l One 
finds that there are in fact two travelling branches of so- 
lutions with the same kinetic energy and opposite phase 
speed. These solutions merge at e = 0.04547 and trans- 
form into a family of nonlinear steady solutions which 
corresponds to a cut at fixed Re of the surface of solu- 
tions in the Re — e — amplitude space. We conjecture that 
this behaviour is a reminiscence in the nonlinear regime 
of the coUision between the first two real eigenvalues of 
the linear problem and their subsequent splitting into 
a complex conjugate pair (see MH05 for full details re- 
garding the Hnear problem). Note finally that Fig. [5] is 
in qualitative agreemenlji^ with Fig. 6 of MHOS. 

The behaviour of steady solutions at large Re suggests 
that solutions with smaller and smaller values of e can 
be found provided that Re is increased sufficiently. How- 
ever, a careful examination of Fig. |4] reveals the presence 
of intriguing "double" bifurcations at large Re (see for 
instance the closed loop in dotted line emanating from 
the filled black circle at Re = 868.2 and e = 0.04547) 
that are not predicted by standard linear stability anal- 
ysis for infinitesimal velocity perturbations with a single 
wavenumber a (this also occurs in Fig. 12 of MH05). 
This behaviour can be explained as follows. Consider the 
equation for the perturbed mean fiow 

Re {dtu + dyuv) — dyU + eRew, (9) 

supplemented by the set of boundary conditions ^ on 
u. Eigenmodes of the linearized version of this equation 
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FIG. 7: Nonlinear continuation of steady a = 0.2 two- 
dimensional solutions of the £-perturbed plane Couette flow 
without the artifact mean flow perturbation as a function 
of e, for Re = 1000 (full line), Re = 2500 (dotted line), 
Re = 9 000 (dashed line) and Re = 30 000 (dash-dotted line). 
(iV,,iV,) = (32,48). 

take on the form 

u = Asm{m:y), (10) 

and their growth rate 7 is given by 

7 = e - (nTr) VRc . (11) 

These Hnear modes have no self-interactions, therefore 
they are also exact steady solutions of the fully nonhnear 
equations for any amplitude A when eRc — (riTr)^. For 
n = 2, we precisely obtain the bifurcation points observed 
in Fig. H at Re = 868.2 for e = 0.04547 and at Re = 
789.6 for e = 0.05, which hints that these mean flow 
modes participate in the dynamics of high Re solutions 
described earlier. Note that the original idea of MHOS 
was to introduce an e perturbation in order to cause the 
least stable eigenmode of Plane Couette flow with = a 
to bifurcate, not to influence the whole dynamics of the 
flow including that of the kx = Fourier mode. 

It is therefore natural to ask whether or not high Re, 
low e branches survive when the artiflcial mean flow 
eigenmode is removed, and if solutions with e = can 
be found subsequently. Thanks to the way the problem 
has been set-up numerically in this study, it is straight- 
forward to get rid of this model artefact by adding the 
e perturbation in Eq. ^ but not in Eq. ([5]). Fig. [6] 
shows that branches of solutions at small Rc computed 
with the corrected model look similar to those in Fig. lU 
which were obtained with the original MH05 model. The 
main difference is for the lowest Re branch represented in 
these figures, which is supercritical in the MHOS model 
whereas it is subcritical when the mean fiow perturba- 
tion is removed. A connection with high Rc branches is 



also found in the corrected model, but in that case the 
branches do not bifurcate at some e-dependent critical 
Re, even though they have a rather low kinetic energy at 
large Re. In order to discover nonlinear plane Couette 
flow solutions, we flnally attempted to continue solutions 
obtained at a given Rc with respect to e. This contin- 
uation is presented in Fig. [7] for different values of Re. 
The amplitude of the solutions is seen to asymptote to 
inflnity in the limit of vanishing e for all selected values 
of Re. The natural conclusion is that no two-dimensional 
nonlinear solutions can be found for plane Couette flow 
using this homotopy technique. 



V. CONCLUSIONS 

The results presented in this study tend to show that 
one cannot identify steady nonlinear two-dimensional so- 
lutions in plane Couette flow using either the Poiseuille- 
Couette homotopy of CEOS or the e-perturbation homo- 
topy suggested by MHOS. We have been unable to recover 
the solutions reported by CE95 and have pointed out 
that they may be spurious. It is therefore also likely that 
the results of Cherhabili and Ehrenstein— suffer from 
the same problem. Other more general homotopies like 
Eq. ([HI) may however be helpful to discover a connec- 
tion between Poiseuille and Couette flow. The steady 
solutions obtained by MHOS at large Rc and small e 
values were shown to be influenced by an artifact of 
the perturbed equations, however using a revised ver- 
sion of their model did not help discovering new solu- 
tions for plane Couette flow because the corrected solu- 
tions found at small e asymptote to inflnite amplitude in 
the £ = limit. This leaves the issue of the existence 
of two-dimensional nonlinear solutions in plane Couette 
flow unsettled. Another interesting avenue of investiga- 
tion would be to consider the fate of the inviscid solutions 
discovered by Kidai^^ in a uniform unbounded shear flow 
when a small amount of viscosity is added. The relevance 
of this type of solutions to wall-bounded flows is however 
quite unclear. 

To conclude, it is worth pointing out that discover- 
ing nonlinear two-dimensional solutions in plane Couette 
flow is also relevant to the problem of the nonlinear sta- 
bility of linearly stable swirling flows. Indeed, plane Cou- 
ette flow with constant (spanwise) rotation can be used 
as a local model of differentially rotating flows. Since 
two-dimensional incompressible dynamics in plane Cou- 
ette flow is not affected by the presence of a Coriolis 
force associated with a rotation vector perpendicular to 
the (x, y) plane, possible two-dimensional nonlinear so- 
lutions in non-rotating plane Couette flow would also be 
solutions in plane Couette flow with constant spanwise 
rotation for any rotation regime, including centrifugally 
stable anticyclonic regimes for which no purely hydrody- 
namic solutions are currently known at alU^. 
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